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We apply analytical and numerical methods to study the linear stability of stripe patterns in 
two generalizations of the two-dimensional Swift-Hohenberg equation that include coupling to a 
mean flow. A projection operator is included in our models to allow exact stripe solutions. In the 
generalized models, stripes become unstable to the skew-varicose, oscillatory skew- varicose and cross- 
roll instabilities, in addition to the usual Eckhaus and zigzag instabilities. We analytically derive 
stability boundaries for the skew-varicose instability in various cases, including several asymptotic 
limits. We also use numerical techniques to determine eigenvalues and hence stability boundaries of 
other instabilities. We extend our analysis to both stress-free and no-slip boundary conditions and 
we note a cross over from the behaviour characteristic of no-slip to that of stress-free boundaries as 
the coupling to the mean flow increases or as the Prandtl number decreases. 

Close to the critical value of the bifurcation parameter, the skew varicose instability has the same 
curvature as the Eckhaus instability provided the coupling to the mean flow is greater than a critical 
value. The region of stable stripes is completely eliminated by the cross-roll instability for large 
coupling to the mean flow. 



I. INTRODUCTION 

Spiral defect chaos (SDC) was discovered in low- 
viscosity convection nearly 20 years ago [l|, and yet 
much of the detail of its origin remain unexplained. The 
Swift-Hohenberg equation (SHE), originally proposed as 
a model equation for fluctuations near the onset of con- 
vection 0] , has yielded considerable insight into the ques- 
tion of wavenumber selection in pattern formation [3|. 
However the SHE has a major drawback as a model of 
low-viscosity convection: mean flows are excluded. The 
importance of the large scale mean flow on the stability of 
convection rolls was first investigated by Siggia and Zip- 
pelius y. Mean flows are known to play an important 
role in the dynamics of spiral defect chaos ^ , and so two 
related generalizations of the SHE have been developed 
to include the effects of mean flows [a, 01 • 

SDC was found numerically in solutions of one of these 
generalized Swift-Hohenberg equations HJ^I and of the 
Boussinesq equations for convection 0, [9l-[l2| . The gen- 
eralized SHE can reproduce some qualitative features of 
spiral defect chaos, at least as transients, though direct 
comparison between the model and convection appears 
not to be appropriate |13| . 

Including mean flows in these models allows an inter- 
esting long-wavelength instability, the skew-varicose in- 
stability (SVI) [ij] whose spatial dependence is neither 
entirely longitudinal nor transverse to the local pattern 
wavevector. The SVI resembles the Eckhaus instability 
but the most unstable modes are those at an angle to 
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the original roll axes. In the SVI, rolls bend and become 
irregular in order to decrease their effective wavenum- 
ber, and often dislocation pairs will form [3. Indeed, 
chaotic spiral patterns have been observed during the 
transition from the conducting state to rolls. Mean flows 
are also crucial to the formation of SDC: Chiam et al. [g| 
showed that spiral defect chaos collapses to a stationary 
pattern when the mean flow is quenched. The onset of 
SDC and defect chaos has therefore been tentatively as- 
sociated with the occurrence of the SVI [l3, [l^ [lg| . It is 
this connection that motivates this detailed investigation 
into the SVI. 

The classical problem of linear stability of convection 
rolls with stress-free horizontal boundaries near the on- 
set of convection has been studied by a number of au- 
thors J17h20| | , while the linear stability of convection with 
no-slip boundaries [21[ is much less well understood. An 
additional problem that occurs in the analysis of the SVI 
is that there is no consistent relative scaling of lengths 
parallel and perpendicular to the roll axes, owing to the 
singular nature of the slow length scale expansion of the 
stability problem, as detailed below. 

Zippelius and Siggia [l3l were the flrst to study the 
SVI. They used asymptotic scalings to derive modula- 
tion equations from the original equations of convection 
with stress-free boundaries, and showed that mean flows 
have an important role. The stability analysis includes 
some restrictive assumptions about the wavenumbers of 
perturbations, and thus does not capture all mechanisms 
of instability [22[ . Bernoff [1^ extended the work of Zip- 
pelius and Siggia by deriving a set of equations for the 
amplitudes of rolls and the mean flow, retaining terms 
with mixed asymptotic orders. The effect of mean flow 
on the SVI for stress-free and no-slip convection has also 
been considered by Ponty et al. j23| . who derived the 



SVI boundary for stress-free and no-slip convection from 
the Boussinesq equations for convection. However, the 
perturbation terms that they included in their calcula- 
tions do not adequately capture the mechanism of the 
SVI in all circumstances [2J]. Milke [2^ carried out a 
more comprehensive study of stability of rolls by means 
of a Lyapunov-Schmidt reduction. He recovered domains 
in Rayleigh, Prandtl and wave number space where con- 
vection rolls are unstable, the so-called Busse balloon. 

In this paper, we revisit the two generalized SH models 
that include a coupling to the mean flow [fl 0, [2a, [23] , 
and analyse in detail how the skew-varicose instability 
depends on the parameters in the models. 

Although the model equations of interest have been 
well studied, a complete analysis of the SVI in the models 
is not available. One difficulty in the analysis comes from 
the contribution of terms proportional to fc^P/(fc^ + P), 
where (fc, /) is the small wavevector of the perturbations 
associated with the SVI. Terms like this are responsible 
for the absence of a consistent asymptotic sca ling in the 
limit of small amplitude and small k and I |28| . This 
difficulty is resolved in this paper. 

The remaining part of this paper is organized as fol- 
lows: in section (II), we present the two generalizations 
of the Swift-Hohenberg equation, both of which incorpo- 
rate the effect of a mean flow. We discuss two operators: 
Va is a projection operator included in the model to allow 
exact stripe solutions, and T^ is a filtering operator that 
suppresses the cross-roll instabilit y. J -'^ is present in an 
original formulation of the models |2g ; Va was suggested 
by Ian Melbourne |20|. In section (HI), we present a de- 
tailed analysis of the linear stability of the stripe solution, 
showing how the SVI can be located analytically. We ex- 
plore the importance of the mean flow in long- wavelength 
instabilities. The structure of the maximum eigenvalue 
for instability for small k and / is also presented. 

In section (IV), we analyze the zigzag and Eckhaus 

long-wavelength instabilities. Section (V) describes the 

interesting skew-varicose instability. We also consider 
various asymptotic limits. 

The work is extended in section (VI) to the stress- 
free boundary condition case, for which there are skew- 
varicose and oscillatory skew-varicose instabilities. Sec- 
tion (VII) covers short-wavelength instabilities (the 
cross-roll instability) for stress-free and no-slip bound- 
ary conditions. Some numerical illustrations of stability 
boundaries with growth rates of perturbations are illus- 
trated in section (VIII), where we also discuss some cu- 
rious behavior of the growth rates of perturbations. Fi- 
nally, in section (IX), we show how the region of stable 
stripes is affected by the coupling to the mean flow and 
how it can be completely eliminated if the coupling is 
strong enough, in the stress-free case. We conclude in 
section (X). 



II. PROBLEM DEFINITION 

In this section, we set out the two models we will inves- 
tigate, and discuss basic properties of the models. Both 
models are generalizations of the two-dimensional Swift- 
Hohenberg equation, where a real field 'ip{x,y^t), repre- 
senting the amplitude of convection, couples to a mean 
flow given by a stream function (^(x, y, t) and its vertical 
vorticity uj{x, y, t) ~ — V^^. 

A. Description of Models 

The standard Swift-Hohenberg equation [2| is 



^ = [m-(i + v^)^]^-V'^ 



(1) 



where 'ip{x,y,t) represents the pattern-forming field, and 
/i is the driving parameter (in convection, fj. represents 
the temperature difference between the top and the bot- 
tom layer), taking the value zero at the onset of pattern 
formation. Both models introduce a (U • V) ip term to 
the right-hand side of ([T]), where U(a;, y, t) is a mean flow 
calculated from the stream function C,[x, y, t): 
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The mean flow has vertical vorticity u!{x,y,t) = — V^^. 
The way that vorticity is generated by nonlinear forcing 
from Tp differs in the two models. 

In the first model [3], the vertical vorticity Lij{x,y,t) 
has its own independent dynamics: 



^ + (U.V)^=[m-(1 + V2)2]v--7'„(V'), 



(2) 



il-P.(V^-c^) 
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(3) 



where Pr^ c and gm are parameters. The Prandtl num- 
ber Pr (the ratio between kinematic viscosity and ther- 
mal diffusivity) is effectively a viscosity parameter for the 
mean flow, which plays a much greater role in low Prandtl 
number convection. Indeed, in the limit of large Pr^ the 
vertical vorticity is hardly excited and the dynamics of if) 
becomes purely relaxational [3|, [ij, [30| , reducing model 1 
back to the SHE. The coefficient Qm is a coupling pa- 
rameter that controls the strength of the mean flow ef- 
fects relative to the ordinary Swift-Hohenberg nonlinear 
term —■0'^. The parameter c models the effect of top 
and bottom boundary conditions on the vertical vorticity, 
with c = corresponding to stress-free boundary condi- 
tions (where uj = constant is an allowed solution of the 
linearised vorticity equation), and c 7^ corresponding to 
more realistic no-slip boundary conditions [uj decays to 
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FIG. 1. (Color online) Growth rates as functions of wavenum- 
ber K = jKj. Dashed blue curve: cri for ^ = 0, which peaks 
at K = Kcriticai ~ 1- Solid green curve: CT2, which takes the 
value —Pre' a.t K = and decreases as K increases. 



zero in the absence of nonlinear forcing) . The operators 
Va and T-y are explained in more detail below. 

The second model 0, |2^, |23| has the vertical vorticity 
slaved to the nonlinear driving term: 



'dt 



+ (U • V) V = [a^ - (1 + V2)2] i, - V^ (V/') , (4) 



= -3^^ [V(V2^0 X VtA] -2, 



(5) 



so the vertical vorticity responds instantly to the non- 
linear driving. The coefficient t; is a coupling parameter 
that controls the relative strength of mean flow effects 
compared to the ordinary nonlinearity. 

We use two operators Va and F^ to ease the analysis 
and to ensure that the SVI is not pre-empted by other 
instabilities. They both act as filters in Fourier space; 
the first is a projection: 



Vaie""^-^) = 



;k : 



if |K|<a; 
if |K|>Q. 



By setting a = 2.5, we allow stripes with a single Fourier 
mode with wavenumber close to one to be exact solu- 
tions of the PDEs dlHSl) and (gHS]) ^. This is shown in 
section (|IIB[) below. The second operator, J^^, reduces 
short- wavelength modulations of the mean flow |26j : in 
Fourier space, the operator is defined by 



^^(e^K-) 



g-7lK|^gjK.x_ 



Throughout this work, we set 7 = 2.5. The effect of this 
filtering on short-wavelength instabilities, particularly 
the cross- roll instability, is discussed in section (|VIII) . 



B. Basic Properties of the Models 



For model 1, linearisation yields: 



Pr (V2 - c2) 



UJ 



only the first equation is relevant to model 2. Normal 
mode solutions to these linear equations are given by 

■0 = ^^gTit+iKx ^^j ^ ^ ^^gaat+iKx^ whcrC (Ti and cr2 

are growth rates, K is a wavevector, and Fi and F2 are 
constants. Substituting these into the linearized equa- 
tions gives (Ti = /.J — (1 — K"^) and 02 = —Pr (^K^ + c^) . 
These are shown in Figure [1] for ^ = 0, when the 
trivial solution is marginally stable. The most unstable 
wavenumber is KcriUcai = 1, and for /i > 0, a band of 
wavcnumbcrs close to X = 1 is linearly unstable, signal- 
ing the onset of pattern formation. The vorticity w (in 
model 1) is always linearly damped, unless c ~ 0. 

We define q = K—K critical and thus the trivial solution 
loses stability for any q at ^Existence, where 



l^Existe 



= {l-{l + q)r 



Note that as g -^ 0, fiExistence -^ M^- 

In this paper, we are interested in the stability of the 
nonlinear equilibrium stripe solution of the PDEs. We 
note that 



A.-(1-(1 + <Z)2)2 



with /? 



(6) 



is an exact solution of both models. This can be shown 
by substituting the expressions for ipo and cjq into the 
PDEs: 

= /i^^ (^e^i+1)- + e-»(i+9)-) _ (1 + v2)2 (^^ (e*(i+«)- 

The nonlinear term in the vorticity equation is zero. 
Next, we note that 



(1 + v2)2e±^(l+9)- = (1 - (1 + g)2)2e±*(l+9)- 



and 



■'2.5 



13 L^iT-+<i> + e-ia+<i)Ay 



3/3 y^ 



(« 



1(1+9)2 



^-i(l+q)x 



provided that q is between —0.167 and 1.5. The conse- 
quence is that (iPqjLUq) is an exact solution. In the limit 
of small q, we have f5 ^ [j-L- 4g2)/3, so stripe solutions 
exist when ji > 4g2. 

The advantage of using the projection Va is that it al- 
lows this exact stripe solution of the PDEs [231. The 
alternative would be to consider the limit of small /i 
and /3, but by having an exact solution, which matches 
the asymptotic result we would obtain without the pro- 
jection, we do not have to be concerned with the relative 
sizes of these parameters compared to other small pa- 
rameters that will be introduced below. 

The two models can be related to each other close to 
onset, regardless of the projection and filtering. By scal- 
ing jjL = C (e^) with e ^ 1, introducing a slow time 



scale d/dt — > e^d/dt, and assuming that the wavenum- 
bers that are excited in the vorticity variable lo are of 
order ©(e), the largest term on the left-hand side of 
the vorticity equation ^ in model 1 is Presto. Thus 
model 1 reduces to model 2 in this limit, with the rela- 
tion g = grn/{Prc^). 



III. LINEAR STABILITY OF STRIPES 

The linear stability theory for stripes in the SHE is 
well known [28|: stripes with wavenumber 1 + q exist 
provided /? > 4g^, and they are stable with respect to 
the Ekhaus and zigzag instabilities provided /? > 12q^ 
and q > 0. Once mean flows are included, the zigzag 
instability needs to be modified at finite Prandtl num- 
ber by the presence of mean-flow modes with non-zero 
vertical vorticity [13, 1221 ■ The Eckhaus instability is un- 
changed. We consider the stability of the stripe solutions 
with respect to long-wavelength perturbations, deriving 
three (model 1) or two (model 2) linear ODEs for the 
perturbation amplitudes, and so determine the parame- 
ter regions in which the stripe configuration is stable, as 
well as the boundary of the SVI. 



A. Linearisation 

We proceed by considering perturbations to the basic 
stripe solution. We suppose that the vorticity pertur- 
bation contains wavevectors {k,l) and {—k,—l). These 
interact with the wavevectors (1 -t- q, 0) and (— 1 — q, 0) 
in the stripe solution to give four new wavevectors, and 

I 



so the perturbed solution can be written as ijj = -00 + "0 
and uj = ujq + uj , with wg = and 



ip = A{t)e 



i{l+q+kd)-x 



B{t)e 



i{-l-q+k,l)-x 



+ A(i)e-''(^+'?+''''')-^ + B(t)e'(i+9-'^^-')"', (7) 
J = C(i)e*('=^')-^ + Cifye^-^^-^^-"" . (8) 

We can calculate the stream function Q from w(a;, y, t) = 
— V^C by inverting the Laplacian, and hence obtain the 
mean flow: 



U 



fc2 + Y' 



(c{ty 



.i{kj)-x 



C{t)e 



-i{k,l)--x 



{l.-k) 



We substitute the expression above into the two models 
and linearize (assuming that A, B and C are small) . Ex- 
amining the coefficients of e'(i+'J+*'''')'^ and e'(~i~9+fe.0x 
results in linear ODEs for A and B: 



A = 



M 



l-[{l + q + kY + 1' 



-6/?M 



-3/3B 



^(i + g)V? 



fc2 



B 



{-l-q + kf + f 



P 
-&fi]B 



C, (9) 



The term k^ + P , which appears in the denominator of 
the governing equations for A and B, arises from invert- 
ing the Laplacian when calculating C and hence U in the 
linearisation of the (U • V) V' term. The equation for C 
differs between the two models. In model 1. the coeffi- 
cient of e*('^'')'' yields 



In Fourier space, the effect of the filtering, F^, is to re- 
duce the amplitude of a Fourier component; for the func- 
tion e*(*'''')'', J%y reduces the amplitude by the multiplica- 
tive factor e~^ '■'' +' K 



Equations P^[TO| and ([12]) for model 2 yield: 



Equations ([9HTT|) for model 1 can be succinctly ex- 
pressed as: 



Ml M2 M3 
M2 Mi -M3 
9mM^ 9m Me Mr 



Ml + gMsM^ M2 + gMsMg 
M2 ~ gM^M^ M4 - gM^Me 

= J2 
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= Jl 
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c 




\ c 



Here, we use the abbreviations: 

Mi^fi-[l-{{l + q + k)^ + P)]'' 

M2 - -3/3, 

M3 = -l{l + q)^/{k^ + f), 



6/3, 



C = -Pr (fc2 + /2 + ^2) C + T^gmKl + q)^/p [[- {I + q - kf - l^ + {1 + g)^] B + \^{l + q + kf + l^ - {1 + g)2] A^ . 

(11) 
In model 2, with no intrinsic dynamics for co, there is an algebraic relation between C, A and B: 

C = J-^.g/(l + q) V^ (((1 +q + kf + f-{l + qf)A + (-(1 + q - kf - l^ + (1 + qf)B) . (12) 

I 



AU = M - [1 - ((-1 -q + kf + P)Y - 6/3, 

Ms = e-'''('='+'')/(l + q)y/P[2k{l + q) + k" + P] 
Me = e-^'e^'+'^'^Cl + g)y^[2fc(l + g) - fc^ - f] 
Mj = -Pr{k^ +1^ +c^). 



We note that in the hmit (fc, /) -^ (0, 0), we have M7 « 
—Prc^, so it is not surprising (looking at the bottom hne 
of the 3x3 matrix for model 1) that long-wavelength in- 
stabilities in model 1 will depend only on the combination 



B. Determinants in the limit of small k and I 

The characteristic polynomials (and hence the eigen- 
values, traces and determinants) of each of these Jacobian 
matrices are even in k and /. Bifurcations occur when an 
eigenvalue crosses through zero. The stripe solution is 
stable only if all eigenvalues are negative for all {k,l), 
so we are interested in extreme values of the eigenvalues 
as functions of k and /. It can be readily checked that 
a zero extreme value of the eigenvalue corresponds to a 
zero extreme value of the determinant. Consequently, we 
use the determinants of Ji and J2 to assist our analysis 
of instabilities. The determinant of Ji, Det(Ji), is: 



where all coefficients P^ etc. are functions of /i, q, Pr, gm, c and the filtering e ' C^ +' ), xhe determinant of J2, 
Det{J2), is: 



(2) 

where all coefficients Pj etc. are functions of /z, q, g and 
the filtering e~^ '■'^ +'"•'. The traces of the two Jacobians 
can be written as 

Tr( Ji) = -6^ - Pr c^ + {A - Pr - 12(1 + qf) k^ 
+ (4 _ P^ _ 4(1 + qf^ f -2{k^ + py 

and 
Tr( J2) = -6,3 + (4 - 12(1 + q)^) P + (4 - 4(1 + qf 
-2/35(1 + g)2e-^'('='+'')) P -2{e + e-f , 

where we recall the relationship in ([6]) between /3 and /.j. 
Note that at this point, no approximations or trunca- 
tions have been made in the linear stability problem of 
the stripe solution, by virtue of having an exact solution. 
Our task is now to work out the most unstable eigen- 
values in the limit of small k and I] this is made more 

I 



I 

challenging by the presence of k^ + P in the denominators 
of the determinants above. 

We note that explicit expressions for eigenvalues of 
Ji are not, in general, analytically attainable (though 
the eigenvalues can be calculated numerically). For 
the matrix J2, in the limit {k,l) — > (0,0), the trace 
is Tr{J2) = — 6/3 and the determinant is zero, so, for 
small (/c, I), one eigenvalue will be —6l3+0{k'^+P), which 
is bounded away from zero for a finite-amplitude stripe. 
The other eigenvalue will be close to zero, approximately 
Det{J2)/Tr{J2). Similarly, in the limit (fc, /) -^ (0,0), 
Det{Ji) = 0, so Ji will have an eigenvalue close to zero 
for small {k,l). Since bifurcations occur when an eigen- 
value is equal to zero, this can be detected in both cases 
by considering only the determinants of the two matrices. 
Hopf bifurcations (see section (jVIl) ) require additional 
consideration. We will expand Det{Ji) and Det{J2) in 
powers of k and /, including the filtering e^^''^^ +''^ in 
the expansion. This yields expressions of the form. 



^ ,, , {Ai.2k^ + Bi.2k^P + Ci,2l^) + {Di,2k'' + Ei,2kH^ + Fi.2k^l'' + Gi.2l^)+0{{e+Pf) 
P>et[Jia) = WTP ' ^ ' 



where in model 2, the coefficients are: 

A2 = 12/3(3(1 + qf - 1) - l&q\l + qf{2 + qf, 



B2 = -24/3(1 - 2(1 + qf) - 1692(1 + qf{2 + qf 

+ 4(1 + qfp (3/3 - Aq{l + qf{2 + q)) g, 
C2 = l2p{Pg{l + qf+q{2 + q)), 



i?2 = 6/3 + 4 + 4(1 + qf ((1 + qf + 2) , 
E2 = 9 [472(1 + qfp (4(1 + qf - 3/3 - 4(1 + qf) 
- 4/3(1 + g)2 ((1 + qf + 1)] + 12 + 8(1 + qf 
+ 18/3 - 4(1 + g)^ 

and in model 1, the coefficients are: 



Ai 



-Prc^A2, 



Bi = -Prc'B2, 
Ci = -Prc'^C2, 
Di = -Prc^D2-PrA2, 

El = -Pr c^E2 + /3(-48 - 84^^ - I68q) + i2q^ 
{q + lf{q + 2)\ 

We will hence refer to the determinants of two matrices 
in the general form with no subscripts: 

Det{J) = [{Ak^ + Bk^l^ + Cl^) + {Dk^ + EkH"^ + FkH^ 

(14) 

The values of F and G arc not needed subsequently. 
Stripes are stable if all eigenvalues are less than zero, 
corresponding to Det{Ji) < and £'et( J2) > 0. To sim- 
plify the presentation, we focus on Det(Ji) since the sign 
of Det{Ji) coincides with the sign of the most unstable 
eigenvalue. 



IV. LONG- WAVELENGTH INSTABILITIES: 
ZIGZAG AND ECKHAUS 

Bifurcation points correspond to parameter values for 
which an eigenvalue has zero real part. We first investi- 
gate how the determinants depend on k and /, and then 
use this information to explore how the bifurcation lines 
depend on the other parameters fj,, q and either g or g™, 
Pr and c. In this section we examine the well known Eck- 
haus and zigzag instabilities, which correspond to pertur- 
bations with I = and fc = respectively. 

In the Eckhaus case, with I = 0, (fT4|) yields Det{J) = 
Ak^+Dk'^ + . . .. For small k, this is positive when A > 0. 
Thus for model 1, instability corresponds to Ai > and 
for model 2, instability corresponds to A2 < 0. 

Accordingly, in both cases, the Eckhaus instability, 
which is independent of the mean flow, occurs when 
A = 0: 



^'Eck 



q^{7q^ + 42q^ + 90q^ + 80q + 24) 



(3^2 + 69 + 2) 



Note that in the limit of (7 — > 0, fisck 

l^Eck ^ '^l^Existence- 



12g2 and hence 



0.075 



0.05 



0.025 




FIG. 2. (Color online) Location of the Eckhaus and zigzag sta- 
bility boundaries in the (q, y,) plane (g < 0), for g — 0.5, 5 and 
50. The Eckhaus boundary derived from j4 = 0, is denoted 
in green thick curve. The dashed red curve is the existence 
boundary. The zigzag boundary (blue thin curves) crosses the 
Eckhaus boundary at g ~ —0.015 when g — 50. The behavior 
for small q is approximately /i = —6q/g. Stripes are stable to 
the right of the Eckhaus and zigzag boundaries. 



when C > 0. Thus for model 1. instability corresponds 
to Ci > and for model 2, instability corresponds to 
C2 < 0. Accordingly, the zigzag instability in both mod- 
els occurs when C ~ 0: 



^zigzag 



^Existence 



3g(2 + q) 
5(1 + 9)2' 



Similarly, in the zigzag case, with A: = 0, ((T^ yields 
Det{J) = CP + Gl'^ -|- . . . , which is positive for small I 



where for model 1 we have identified g = g„i/ Prc^ . 

Unlike the Eckhaus instability, the zigzag instability is 
affected by the mean fiow. Vorticity and mean flows act 
as a stabilizing influence on the zigzag instability, which 
is suppressed for larger values of g, resulting in a larger 
region of stable stripes for q < in the (/i, q) stabil- 
ity diagram. Figure [2] shows how the zigzag instability 
boundary behaves for different values of g. Note that 
for large enough (7, it no longer forms the lower stability 
boundary except for very small /i. 

The zigzag and Eckhaus instabilities cross in parame- 
ter space when g and q are related by 



_ 3(3(1 + 9)2-1) 
^ 4g(l + q)4(2 + g)' 

where q < Q. Furtherinore, when g — > 0, we recover 
the standard result for the SHE without mean flow, and 
when 5 -> cxj, we have Hzigzag -^ fJ-Existence- However, 
even in this limit, for small ^ and q, fj-zigzag is approxi- 
mately —Qq/g, so the zigzag stability boundary emerges 
from {q,fi) = (0,0) in a straight line of slope —6/g, as 
can be seen in flgure [2] 



V. LONG- WAVELENGTH INSTABILITIES: 
SKEW- VARICOSE 

The skew- varicose instability, which is driven by the in- 
clusion of mean flow, the strength of which is determined 
by 5, is associated with modes for which the maximum 
positive growth rate occurs when fc 7^ and / ^ 0. Two 
conditions arc required to characterize the SVI: the de- 
terminant should be zero and should have maximum or 
minimum value (for model 1 and model 2 respectively) 
for fc 7^ and I ^ 0. We first express these conditions for 
the SVI in terms of the coefficients A — G oi the expres- 
sion ()14|) . the power scries expansion of the determinants 
of Ji and J2, which we denote simply by Det. We then 
express these conditions in terms of the parameters /x, q 
and either g or gm, c and Pr, in order to locate the SVI 
boundary in the (/^, q) plane. 

There are two different manifestations of the skew- 
varicose instability. In the first (case I), the instability 
emerges from fc = Z = 0, as illustrated in figure [S] In this 
case, the determinant is negative for {k,l) close to (0,0), 
corresponding to A < and C < 0. If we suppose in the 
first instance that we can write the leading order terms 

m the determinant as: JJet ~ PTP — ^ ^ 

then imposing Det = 0, along with dDet/dk = and 
dDet/dl = 0, results in (B^ - AAC){A - B + C) =0 
and p- = -gEf^- In order for k"^ /P to be positive, we 
need B > max{2A, 2C), which excludes B = A + C and 
B = -VIAC. We conclude that 
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A 
(15) 
is the condition for the SVI. However, the truncation 
of Det above is degenerate: the conditions arc satisfied 
along a line in the {k, I) plane, rather than at a point. 
This degeneracy is resolved by restoring the higher order 
terms, as illustrated in figure [3] 

The other possibility for the skew- varicose instability 
is that it can accompany the Eckhaus instability. This 
occurs when A > and C < 0: for I = and D < 0, 
Det is positive for a range of k and attains its maximum 
on the / = axis at a finite k = k^ax- In the {k,l) 
plane, {k^ax, 0) can either be a maximum or a saddle, as 
illustrated in figure 21 We define the SVI (case II) to be 
the point at which (fc,„a2;7 0) changes from a maximum 
to a saddle; at this point the maximum eigenvalue moves 
off the k axis. 

Unlike in the previous case, the growth rate at the 
SVI is positive, since stripes are already Eckhaus unsta- 
ble. Therefore the SVI occurs when there is a degenerate 
maximum at / = and k = kmax, about to become sad- 



dle. The conditions ^M^ = and 



dDet 



at this point 



yield the parameter values at which this variant of the 
skew-varicose instability occurs. 

At I = 0, the first condition implies fc^^j^. = ^^ + 
0{A'^), and the second condition implies {B — A) + {E — 



FIG. 3. (Color online) Case I of the SVI: behavior of the Det 
in the {k,l) plane for A = —0.1 and C = —0.1 {A < and 
C < makes stripes Eckhaus and zigzag stable). The other 
coefScients are D = —1, E = 2, F = —1 and G = —1. (a) 
B = 0.195, giving stable stripes (B^ < 4AC). (b) B = 0.205, 
giving stripes that are unstable to the SVI {B^ > 4AC). The 
positive maximum of the determinant emerges from (fc, I) = 
(0, 0) but occurs with A; 7^ and Z 7^ 0. A negative value of 
the determinant is indicated by black contours while zero and 
positive values of the determinant are in red(gray). 



^) \ 2DJ ~ " "•" "^ ~ '^m.ax 

A, the condition for case II of the SVI is 



at fc 



Hence in the limit of small 



B= (2±E\a + 0{A^), A>OandL'<0, 

with fc2 == ^ + 0{A^) and / = . (16) 

Contours of Det on either side of this boundary are il- 
lustrated in figure m Moreover, the conditions (flS)) and 
(|16l) arc summarized in the schematic diagram in figure [SI 
which also indicates the regions affected by the SVI and 
the point {A, B) = (0, 0), where the two cases coincide. 

We then translate the conditions P^ and (|16p into 
bifurcation lines in the (/i, q) plane. In order to derive 
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FIG. 4. (Color online) Case II of the SVI: behavior of the Det 
in the (fc,Z) plane for A = 0.005 and C = -0.1 (^ > and 
C < makes stripes Eckhaus unstable but stable to zigzags) . 
The other coefficients are D = —1, E — 2, F — —1 and G — 
— 1. (a) B = —0.003, giving SV stable and Eckhaus unstable 
stripes (B < ^^y-Aj ; the maximuni occurs with I = 0. (b) 
B — —0.001, giving stripes that are unstable to both SV and 
Eckhaus instabilities (B > ^^y-A); the maximum moves off 
axis, and {kmax,0) is now a saddle. A negative value of the 
determinant is indicated by black contours while zero and 
positive values of the Determinant are in red (gray) . 



bifurcation lines we have used a branch-following pack- 
age, MATCONT IHj. An example of the calculation for 
model 2 is given in figure [SI which was computed work- 
ing directly with numerically determined eigenvalues and 
computing -t^ and -Mj numerically. In case I, the SVI 
boundaries derived using condition ()15p coincide with the 
numerical computation. However, in case II, condition 
psp agrees with the numerical computation only when 
A is close to zero, as would be expected. The transition 
from case I to case II occurs at {fi,q) = (0.1367,0.2387) 
for 5 = 0.5: for smaller fi, the stability region of stripes 
is bounded on the right by the Eckhaus instability, while 
for larger /i, the Eckhaus instability is preempted by the 
SVI. For g ~ 5, the SVI pre-empts the Eckhaus instabil- 
ity for all fj,. 
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FIG. 5. (Color online) Schematic diagram of the conditions 
for the Eckhaus and SV instabilities in cases I and II in {A, B) 
plane. The Eckhaus instability occurs when ^4 = (dashed 
red line). The blue curve shows the SV stability boundary: 
B'^ = 4AC in case I (A < 0), and B = {^j^A) in case II 
{A > 0). At the intersection of the two stability boundaries, 
{A,B) = (0,0). 
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FIG. 6. (Color online) Numerical computation (in model 
2) of the SV stability boundary in the {fi, q) plane. For 
5 = 5, the SVI pre-empts the Eckhaus instability for all 
H and hence the region of stable stripes is bounded by the 
skew-varicose instability curve. However, for g = 0.5, the 
region of stripe stability is bounded by the skew-varicose in- 
stability curve only when /.t > 0.2387. The crossing point, 
(0.1367,0.2387), of the two boundaries is denoted as a red 
square. For < /x < 0.2387, the Eckhaus precedes the SV 
curve, which reaches the origin as a parabola p — 9.33g^. The 
green thick curve denoted by E is for the Eckhaus boundary 
whereas the dashed red curve is the boundary of existence of 
stripes. 



A feature of the SVI is that it makes stripes unstable 
only for g > and, for any given value of /x, g must 
be large enough for the SVI to preempt the Eckhaus in- 
stability. The crossing point between the Eckhaus and 
skew-varicose instabilities for a given value of q occurs 
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when g < QcriUcai ■ A detailed presentation of this asymp- 
totic result will be discussed in section (|V Ap . The dis- 
tinction between g > gcriUcai and g < gcriUcai is iUus- 
trated in figure IHl where we present SVI boundaries for 
g = 5 > gcriUcai and g = 0.5 < g critical- 

Interestingly, for a fixed g, 5 — > implies ^sv ^^ 
^- Existence- Therefore when g — > 0, the SVI boundary 
coincides with the existence curve of stable stripes for 
small ^. 

The Eckhaus and SVI boundaries are often very close, 
so we present our result in an alternative way in figure [T] 
Figures 7(a) and |7(b)| present the SVI boundaries for 
5 = 5 and g = 0.5 in the [ii/ ii Existence , \ogif){q)) plane. 
This is a better way of illustrating the regions of stable 
stripes (shaded regions) and the behavior of the Eckhaus 
and SVI boundaries as g — ?► 0. Moreover, it shows how 
the coordinates axes A and B from (14) can be defined 
near the SVI-Eckhaus crossing point. 



A. Asymptotic analysis of the SVI boundary 

In this section we focus on the asymptotic behavior of 
the SVI boundary in model 2, in the two cases discussed 
above, g > g critical and g < g critical- Let us consider 
the case when g > gcriticai, where the SVI boundary 
pre-empts the Eckhaus instability boundary. We use the 
condition B^ — AAC = 0, which can be written as 



Fig^ + F2g + F3 = 0, 



(17) 



FIG. 7. (Color online) Numerical computation (in model 
2) of the SVI boundary as a function of fi/fiExistence and 
logio(g). fiExistence = (1 - (H- q)'^Y - The Eckhaus bound- 
ary, yl = 0, is denoted by a green thick curve and the shaded 
region corresponds to stable stripes, (a) p = 5 > Qcriticai- 
Here, [s^sv I P^Exiatance. — >■ 3 as (J — > 0. The SVI precedes 
the Eckhaus for all q values, (b) 3 = 0.5 < gcruicai- Here, 
fJ-sv / f^Existence -^ 2.3333 as Q — >■ 0, which in turn becomes 
fj,sv — >■ 9.3333g^ as g — ^ 0. The point of intersection of 
the SVI boundary with the Eckhaus boundary is denoted by 
{A, B) = (0, 0). A schematic illustration of the A and B axes 
at the crossing point is also shown (see figure [5]). 



for some g, say gEck, which is given by 



QEck ~ n 



3(l + g)2-l 



{l + qf 



The condition ^ = could be used to express gsck as 
a function of ^ if desired. In the limit /x — > and 
g ^ 0, gEck goes to 0.75, which we call gcriticai- The 
expression for gEck and the value gcriticai = 0.75 is 
the same in models 1 and 2 provided we use the rela- 
tion g = g,n/{Prc^). The Eckhaus instability precedes 
the SVI for some range of /x only if g < gcriticai- We 
have found that the SVI boundary for g > gcriticai ap- 
proaches the origin as /i = 12q^, as does the Eckhaus 
curve, whilst it approaches as /x = ng^, with 4 < n < 12 



where the i^^'s are functions of fi and q. In the limit of 
very small /i and g, 

Fi « i(16384g'^ - 8192^g'^ + 1024^*^^2 _ 256^*^? + 16/), 

(18a) 

(18b) 

(18c) 



F2 « -(-24576g^ + 8192^iq^ - 512^i^q - 64//^), (18b) 
F3 w (9216g'' - 1536Mg^ + 64/1^). 



We note at this point that equation p?)) is valid for 
model 1 (with the same values of Fi, F2 and F3) when 
we identify g with gm/Prc^. For smallish g > gcriticai, 



as /z — > 0, /x '--^ g as shown in figure 7(a) and hence 



functions in equation p8|) can be taken as Fi ^ q° , F2 ^ 
q^ and F3 ^ q^ . Therefore the equation p?|) can be 
approximated as F^ = 0, in which case /i = 12g^. This 
can be improved by including F2 and hence i*25 + ^3 = 
in which case we obtain 



^^sv 






9 - 24gg; 



{g > gcriticai , q < i/g) (19) 



and hence /i — > 12g^ when g — > 0. Therefore when 
g > gcriticai 1 the SVI boundary for small fj, has the same 
curvature as the Eckhaus boundary. Note that figure 7(a) 
shows how the SVI boundary for g = 5 coincides with the 
Eckhaus boundary as g — >■ 0. The limit of very large g 
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FIG. 8. (Color online) Numerical computation of the SVI 
boundary on a logarithmic scale for g = 10' for i — 1,2, 3, 4, 5 
and 7. The transition from /i ~ g to /i ~ q'^ occurs when 
g ex 1/g. In the limit of small q, the SVI curve is tangent 
to ^ = I2q^ , whereas in the limit of large g, the SVI curve 
goes as /i ~ 8g as /i increases. Both asymptotes, /i ~ 8g and 
H ~ 125^ are denoted by red thick lines. 



fJ-SV 



Qcritical 



g+l\q , {g < QcrtUcal , 9 < 1). 



(21 



This relation is shown numerically using MATCONT [31| 
to be a very good approximation; fJ-sv / fJ-Existence be- 
haves linearly with g with a gradient 8/3. In addition, 
the numerical simulation illustrated in figure mbjl shows 
at g = 0.06, ^/ ^Existence ~^ 2.6 as g — > which agrees 
well with equation ((2T|) . All these explicit results arc for 
model 2. When g > gcriUcah the expressions for con- 
ditions given in (|15p are the same for model 1 with the 
relation g = gm/Prc^. Therefore, equations p^ and 
pO)) are the same for model 1. When g < gcriUcai, for 
model 1, the condition for the SVI given by equation (|16p 
with g = gm/Prc^ is not the same as in model 2. How- 
ever, for fixed Pr and c 7^ 0, we find that fJ-sv/fJ-Existence 
behaves linearly with g^ and hence equation (j2ip holds 
for model 1. The case where c = is considered in the 
next section. 



is also of interest, as it corresponds to stress-free bound- 
ary conditions (see section (jVI[) below). In this limit, we 
expect the SVI boundary (as shown in figures [14] and [TBI 
below) should have fi ^ q [l9[. From equation (1181) with 



fi ^ q, wc have Fi ^ q'^, F2 ^ q'^ and F3 ^ q^ and 
Fig^ + F25 + F3 - g2 (^2^2 ^ q^ _^^y If 95 < 1, we 
set F3 = and recover n = 12q^. li qg ^ 1, wc set 
-Fi = and obtain ^ ~ 8q. Again, this can be improved 
by setting Fig + F2 = 0, and we obtain 



fJ'SV 



^qg + 3 
2qg~3 



, {g > gcriUcai , q > 1/5)- (20) 



Finally, we note that the transition between /i ^ 12q^ 
and /i ~ Sq will occur when qg is of order unity, so 
Qtransition ^ -■ Thcsc three regimes are illustrated in 
figure m 

In the skew- varicose mechanism for g > gcriUcai, the 
maximum of Det is attained for perturbations of a mode. 
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and L, 



~ \/^ and so 



— )• 0. As shown 
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Secondly, let us consider the case g < gcriUcai, for 
which the SVI boundary lies in the Eckhaus band for 
small fi. We do not have an exact criterion for the SVI 
in this case, though equation (16) is an approximate cri- 
terion. However, for small /i, we expect fisv to depend 
approximately linearly on g. We know from equation 
(19) that for g = gcriticah IJ-sv ^ 12g^ as g -^ and 
from conditions in equation (16) in the limit g — > 0, we 
have ^sv — > 4g^ as g — >■ 0. A linear interpolation bc- 



VI. LONG-WAVELENGTH INSTABILITIES: 

SKEW- VARICOSE AND OSCILLATORY 

SKEW- VARICOSE WITH STRESS-FREE 

BOUNDARY CONDITIONS. 

We now consider the case that models convection with 
stress-free boundary condition. In model 1, c is the pa- 
rameter that accounts for the boundary conditions at the 
top and bottom, and stress-free boundary conditions cor- 
respond to c = 0. Using the relation g = g^jPrc? to 
connect the two models, taking the limit 5 — > 00 corre- 
sponds to stress- free boundary condition in model 2. 



When 



0, for any coupling constant ^m, the SVI 



boundary always pre-empts the Eckhaus boundary in 
the (/i, g) plane. To show this we can focus on model 
1, for which the SVI condition is F\g^ ■\- F^gmPT (? 4- 
F^(Pr(?Y ~ with Fi, F2 and F3 given in equation 
P^ : this confirms the relation g = g„i/Prc^. Hence 
c = implies Fi ~ 0, and this gives the criterion for 
instability as /x = 8g for any g^, as above. A remarkable 
property of the SVI in model 1 with stress-free boundary 
conditions is that the stability boundary is independent 
of Pr and gm ■ 

Another instability of interest in the stress-free case 
is the oscillatory skew-varicose (OSV) instability, which 
consists of a long-wavelength transverse oscillations of 
the stripes that propagate along their axis ]1^. For OSV 
modes, the associated eigenvalues are complex. This in- 
stability docs not occur in model 2 because Tr(J2) < 
for small k and I and so complex eigenvalues arc not pos- 
sible. In model 1, the OSV instability does not appear 
for c^ = 2. With c = 0, the regions in the (fc, /) plane 
with positive growth rate and nonzero frequency emerge 
from (fc, I) = (0, 0) in a way that resembles the contours 
for the SVI, shown in figure [31 
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FIG. 9. (Color online) The location of the OSV instability 
boundary for model 1 for c = (stress-free boundary condi- 
tions), Pr = 1 and g^ = 1000, 100,25,5 and 1. Stripes are 
OSV unstable to the left of the OSV boundary. For small /i, 

the boundary is asymptotic to /i = ( ^'^t ^^ ) igm- The Eck- 
haus boundary is denoted in green thick and the existence 
curve is in dashed red. 



FIG. 10. (Color online) Stability diagram in the neighborhood 
of ^ = for model 1 with c^ = 2 (no-slip boundary condi- 
tions), Pr = 1, gm = 50 and 7 = 2.5. Stable stripes are in 
the region indicated, bounded by Eckhaus instability (green 
thick curve) from below and by the SVI boundary (blue thin 
curve) from above. 



We derive the asymptotic behavior of the OSV insta- 
bility boundary. At the point of instability, the eigen- 
values are purely imaginary; this gives two conditions, 
C — AB = and i? > 0, where the characteristic equa- 
tion of Ji is A^ + A\^ + B\ + C = Q. The coefficients A, 
B and C are functions of g™, q, fi, k and I. We have set 
Pr = 1 for illustrating this calculation. For small k and 
I, the condition C — AB = gives. 



Gi 



9m q 



Go 



9mq 



G3 - 0, 



(22) 



after maximizing C — AB over (fc, I). The G^'s are func- 
tions of /i and q. In the limit of very small fi and q, we 
find. 



G,.i 1024 -8192 (^ 



16384 




(23a) 



512- ^r^U 46421 r^^ - 76459 r^ 



G2 



G^w 256 -8192 1 — 



(23b) 
90112( — j -393216[ — 

589824 ( — ) , (23c) 



where we have dropped terms that can be shown to be 
smaller than those retained. 

We note at this point that equations ([22]) and ([23]) with 

2 



512(^^ 



256 = 0, 



(<zVm) « 1 give ^ [^ 

which implies fi — I "'^^ ^ J qgm- Hence, the OSV insta- 
bility boundary has a linear relationship between q and fi 



for small /i, and it bounds the region of stable wavenum- 
bers for negative q. Stripes are stable on the right of the 
boundary. In this asymptotic limit, the point of max- 
imum growth rate in the {k,l) plane can be found at 
the point of maximum of G — AB; this point satisfies 
l/k ~ \/5- Figure M shows the behavior of the OSV in- 
stability boundary in model 1 with c ~ 0, for different 
values of the coupling constant gm. 

When c is increased from zero, the OSV instability 
turns into the so-called oscillatory instability jS^]. This 
has the nature of an oscillatory cross-roll instability, set- 
ting in with non zero k and /. The boundary of this 
oscillatory instability emerges from /3 = 0, the existence 
curve. However, this instability is prominent only for 
Pr c^ ^ 1; for higher values of Prc^, the instability 
moves to larger negative q. 



VII. SHORT- WAVELENGTH INSTABILITIES: 
THE CROSS-ROLL INSTABILITY. 

The cross-roll (CR) instability is so-called because the 
fastest growing disturbances appear to take the form of 
stripes perpendicular to the basic steady stripe pattern: 
these disturbances have non-zero k and / in the limit fi — >■ 
[301 . In contrast to the oscillatory cross- roll instability 
discussed above, the most unstable eigenvalue at the CR 
instability is zero. 

We will show numerically in section (|VIII|) below that 
the CR instability only forms a boundary of the region 
of stability of stripes if g is large. The filtering J'^, dis- 
cussed in section ([H]) , was introduced to suppress the CR 
instability in favour of the SVI, so the location of the CR 
instability boundary depends on 7, while the SVI is not 
influenced by 7 for small enough ^ and q. Since the CR 
instability sets in with non-zero k and I even with small /i 
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FIG. 11. (Color online) Stability diagram in the neighborhood 
of /i = for model 1 with c^ = 2 (no-slip boundary condi- 
tions), Pr = 1, Qm ~ 1000 and 7 = 2.5. Stable stripes are 
in the region indicated, bounded by the SV and CR instabil- 
ities from above, and by the zigzag and Eckhaus instabilities 
from below. Stable stripes exist at g = for range of /i close 
to zero. Growth rates as a function of k and I at the points 
indicated by (a)-(f) are given in figure [T2] 



and q, asymptotic analysis of the type carried out above 
cannot be done. Numerical examples are shown in the 
next section. 



VIII. NUMERICAL ILLUSTRATION OF 
COMBINED STABILITY BOUNDARY. 

We now present numerical results obtained with MAT- 
CONT 



31 



Details of the specification of the conditions 
on the eigenvalues for each of the six types of instability 
are given in the Appendix. We choose two illustrative 
parameter values: gm = 50 and gm = 1000, and we con- 
sider no-slip (c^ = 2) and stress-free (c = 0) boundary 
conditions. 

Figure [TU] shows the stability diagram for model 1 with 
parameter value c^ — 2 (corresponding to no-slip bound- 
ary conditions), Pr ~ 1, g^n = 50 and 7 = 2.5. The re- 
gion of stable stripes is bounded by the SVI from above 
and the Eckhaus instability from below. The zigzag in- 
stability boundary for this value of gm hes below the left 
Eckhaus boundary except for very small /i and q (as in 
figure 2). The CR instability does not occur in this range 
of parameters. 

Figure [TT] and [T51 show the stability diagram for model 
1 with parameters c^ = 2 (no-slip), Pr = I, gm ~ 1000 
and 7 = 2.5. Close to fi ~ (figure [TT|) . a neighborhood 
of (7 = is in the stable regime and the zigzag and the 
SV instabilities bound the region of stable stripes. For 
/i > 5 X 10~^, the region of stable stripes is bounded by 
the CR instability from above and by the Eckhaus insta- 
bility from below. For this value of <?,„ , the CR instability 
boundary crosses the SVI boundary, and the zigzag in- 
stability boundary has a linear behavior for small /i. 

For the same parameter values, the stability diagram 



for a larger range of q and fj, is shown in figure [131 The 
region of stable stripes is bounded by the Eckhaus in- 
stability from below and the CR instability from above. 
The zigzag instability boundary lies close to the existence 
curve and is of less interest for this value of gm ■ 

Figure [T^] shows the eigenvalue behavior at selected 
points in the (q, fi) space from figure 11, as a function of 
(/c, /), showing how stripes can be unstable to one or both 
of the SV and CR instabilities. We note that the CR in- 
stability occurs for reasonably large values of fc « 0.04 
and / « 0.2. In contrast, in the SVI, contours of posi- 
tive growth rate emerge from {k,l) = (0,0). When both 
instabilities exist, two separate peaks of growth rates ap- 
pear. For large g, these contours can join to form one 
large contour. 

We have computed the stability diagrams for model 2 
with g = 500 and g = 25, and these are qualitatively 
the same as [10] and [11] consistent with the relation g = 
gm/Prc^. 

Figure UM and \W\ similarly show the instability bound- 
aries for model 1, with parameters c ~ (corresponding 
to stress-free boundary conditions), Pr = 1, 7 = 2.5 and 
gm = 50 and gm = 1000. The results agree qualitatively 
with earlier calculations by Bernoff [l9[, who found a 
similar linear relation between // and q for the SV and 
OSV instabilities in convection with stress-free boundary 
conditions; he did not consider the CR instability. 

Disregarding the CR instability, stripes would be sta- 
ble between the OSV instability and SVI boundaries. 
However, as seen in figure [T31 the SVI is always pre- 
empted by the CR instability; these boundaries appear 
to be parallel for larger /i. For /i < 0.32, there are no sta- 
ble stripes. For higher fi, the stable region is bounded by 
the CR instability from above and by the OSV instability 
from below. 

Figure [TS] shows the change of structure of the eigen- 
values when moving from left to right in the stability 
diagram shown in figure [TJ] At ^ = 0.1, we selected 
four different wavenumbers: q = —0.05, where stripes 
are OSV unstable, q = —0.02, where stripes are CR and 
OSV unstable, q = 0.0058, where stripes arc CR unstable 
but OSV stable and q = 0.05, where stripes are CR and 
skew-varicose unstable, though the distinction between 
these two instabilities has become blurred. 

Figure [THJ presents instability boundaries for stress- free 
boundary conditions with Pr ~ 1, gm = 50 and 7 = 2.5. 
This provides an illustration of the change of the CR 
instability boundary with gm- For gm = 50, the CR 
instability boundary crosses the SVI boundary and the 
effect of the CR instability is reduced. 



A. Curious behavior of growth rates 

In previous sections, we considered the growth rate 
behavior of {k,l) close to zero. However, when a larger 
range of k and / is considered, we noticed additional insta- 
bilities in the regime where stripes are already unstable. 
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FIG. 12. (Color online) Growth rates of perturbations at selected {fJ-,q) indicated by (a)-(f) in figure [TT] (a) /i — 8.5 x 10^'', 
q = 0.002; stripes are CR unstable, but growth rates for k and I close enough to zero are negative, (b) fj, = 5.5 x 10^^, q — 0.003; 
stripes are CR and SV unstable. Growth rates for k and I close enough to zero become positive due to the SVI and there are 
two distinct zero contours, (c) fi = 5.5 x 10"*, q = 0.0034; stripes are CR and skew- varicose unstable. One large contour 
encloses both the SV and CR instabilities, (d) fi — 2.5 x 10"**, q = 0.00225; stripes are skew- varicose unstable, and growth 
rates for k and I close enough to zero are positive for a range of polar angles, (e) fi = 2.5 x 10"*, q = 0.0025; same as case (b), 
but the maximum occurs in the SV region, (f) fj, = 2.5 x 10"'*, q — 0.003; same as case (c), but again the maximum occurs in 
the SV region. The zero contour is denoted in black (outer) and contours of positive growth rate are denoted in red (inner in 
gray). 
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FIG. 13. (Color online) Stability diagram with parameters as 
in figure [11] covering a larger range of q and /i. Stripes are 
stable in the shaded region, bounded by the Eckhaus and CR 
instabilities. 



These instabilities occur even in the SHE and so are not 
related to presence of the mean- flow. They have not been 
studied before, though they are of less interest since they 
do not bound the region of stable stripes. We consider 



FIG. 14. (Color online) Stability diagram for model 1 with 
c = (stress-free boundary conditions), Pr = 1, gm = 1000 
and 7 = 2.5. Stripes are skew-varicose unstable to the right 
of the SVI boundary, OSV unstable to the left of the OSV 
boundary and CR unstable to the right of the CR boundary. 
Stable stripes exist for /i > 0.32; the region of stable stripes 
is shaded, and is bounded by the OSV instability boundary 
on the left and the CR boundary on the right. Growth rates 
as a function of k and I at the points indicated by (a)-(d) are 
given in figure 15. 
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FIG. 16. (Color online) Stability diagram for model 1 with 
c = (stress- free boundary conditions), Pr = 1, gm = 50 and 
7 = 2.5. Stable stripes exist for /.t > 0.02; the region of stable 
stripes is bounded by the OSV boundary on the left and the 
CR boundary and then the SV on the right. 



(c) 



(d) 



FIG. 15. (Color online) Growth rates of perturbations at 
selected (/i, q) indicated in figure [Til all with /j. — 0.1. (a) 
q = —0.05, where stripes are OSV unstable, (b) q — —0.02, 
where stripes are CR and OSV unstable, (c) q = 0.0058, 
where stripes are CR unstable but OSV stable, (d) q — 0.05, 
where stripes are CR and skew-varicose unstable, though 
the distinction between these two instabilities has become 
blurred. The zero contour of the real part of the eigenvalue is 
denoted in black (outer) and contours of positive growth rate 
are in red (inner in gray). The eigenvalues in the OSV case 
are complex. 



them briefly here. 

Contours of the growth rates of perturbations for se- 
lected parameter values are shovifn in figure [TTl As shown 
in figure 17(a)[ for a fixed small /x, when q is increased 
from the left existence boundary (g < 0), modes ap- 
proximately in an annulus of unit radius centered at 
{k, /) = (!-!- q, 0) have positive real eigenvalues in addi- 
tion to unstable modes for small k and I, corresponding 
to the Eckhaus and zigzag instabilities. This annulus dis- 
appears for larger q leaving only the Eckhaus and zigzag 
uns table m odes close to {k,l) = (0,0), as shown in fig- 
This process is in part reversed for q 



17(b) 



> 

going towards the right existence boundary: here modes 
are Eckhaus unstable only, but for large enough q, the 
annulus of unstable modes reappears. The behavior of 
growth rates in the right Eckhaus band is illustrated in 
figure [T7(c)] and [T7(d)| 



We have confirmed that this curious behaviour is not 
a result of using the projection operator Pa in the Swift- 
Hohenbcrg equation. 



IX. ANALYSIS OF THE ROLE OF 
MEAN-FLOWS 



We now illustrate the stability diagrams in a different 
manner. For a fixed fi, we show how the coupling to the 
mean flow affects the region of stable stripes. This choice 
of presentation provides useful information for numerical 
simulations of the PDEs in large domains. Figure [18] 
represents the region of stable stripes for model 2 in the 
(g, q/y/JI) plane for /i = 0.1. We choose q/y/JI as the coor- 
dinate for ease of comparison between different values of 
/i. Figure [Is (a) shows the stability diagram for small g, 
where the region of stable stripes is bounded from above 
by the Eckhaus instability for g < 0.574 and by the SVI 
for g > 0.574, and by the zigzag instability from below. 
Figure [T8(b)| shows how for large g, the region of stable 
stripes is bounded by the CR instability from above and 
by the Eckhaus instability from below, and the region is 
eliminated when g >2 x 10"*. Figure \W\ shows the loca- 
tion of the CR instability for ^ == 0.01, 0.001 and 0.0001. 
The upper bound on g beyond which there are no stable 
stripes initially decreases with fj. and then increases. The 
behavior of model 1 is qualitatively the same. 

Figure [20l presents stability diagrams for model 1 with 
c = (stress-free boundary conditions), for /x = 0.1 and 
/i = 0.01. The SV and the OSV instabilities bound the 
region of stable rolls from above and below and the CR 
instability makes the upper bound on g„i. Stable stripes 
exists only for small gm and the upper bound on gm fur- 
ther reduces with /i. 



X. CONCLUSION 

In this work we considered the consequences of includ- 
ing the mean-flow on the generalized Swift-Hohenbcrg 
models. We analysed two models: in the first model 
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FIG. 17. (Color online) Growth rates of perturbations as a 
function of {k,l) for the Swift-Hohenberg equation ^, with 
ip^ replaced by Pa{ip^)- Figures (a) and (b) correspond to 
points in the left Eckhaus band, q < 0, whereas (c) and (d) 
correspond to points in the right Eckhaus band, g > 0. (a) 
Stripes are zigzag and Eckhaus unstable and the unstable 
modes lie approximately in a annulus of unit radius centered 
at {k,l) — (1 + g, 0) have positive real eigenvalues, (b) Stripes 
are zigzag and Eckhaus unstable, (c) Stripes are Eckhaus 
unstable, (d) Stripes are Eckhaus unstable and in addition, 
an annulus of unit radius centered at {k,l) = (1 + q,0) has 
positive real eigenvalues. Thick black: zero contour, thick 
Red (inner in gray): positive contour. Dotted curve: circle of 
unit radius centered at {k, I) = (1 + q, 0). 



vorticity has its own independent dynamics Q- In the 
second, vorticity is directly slaved to the order param- 
eter [2y|. These two models are related to each other 
through g = gm/{Pr c^), where Qm and g are the cou- 
plings to mean flow in model 1 and model 2 respectively. 
Two boundary conditions were considered in this work: 
stress- free (c = 0) and no-slip (c^ = 2). 

In order to explore long-wavelength instabilities, we 
carried out a complete linear stability analysis of stripes. 
We expressed the relevant determinants as power series 
in k^ and Z^, where (fc, I) is the perturbation wavevector. 
We were able to derive explicit expressions for the largest 
growth rates in most cases. This has led to an improved 
understanding of the instabilities of stripes. Unlike in 
previous work [1^, [23, [2^ , we have not had to make as- 
sumptions on the relation between k, I and the amplitude 
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FIG. 18. (Color online) Stability diagrams in {g, q/^/JI) plane 
for model 2 with /i = 0.1 and 7 = 2.5. (a) For small g: the 
region of stable stripes is mainly bounded by zigzag and SV 
instabilities, (b) For large g, the region of stable stripes is 
bounded by the Eckhaus (thick green) and CR instabilities. 
These instability boundaries cross around g — 2 x 10'*, elimi- 
nating the region of stable stripes. The region of stable stripes 
is hatched. 



of the basic stripe solution. This approach has been made 
possible through the use of the projection operator, Pa, 
which allows the exact stripe solution to be written down 
easily ^. 

The skew- varicose instability, which was our main con- 
cern, has two different behaviors: if stripes are stable to 
the Eckhaus instability, in the limit oi fi = 0, the SVI 
goes as /i ~^ 12?^, provided g > 0.75. The most unstable 
wavevector satisfies k'^/P = 0(1). For g < 0.75, the SVI 
boundary crosses the Eckhaus curve, and in the limit of 
/i = 0, it goes as ju ^ aq^ with 4 < a < 12. In model 1, 
the critical g^ is 0.75Prc^. In the large g limit (that is, 
for very low Pr, or for stress-free boundary conditions), 
there is a transition of the SVI boundary from /i = 12q^ 
to /i = 8(7 at a wavenumber satisfying q oc l/g. 

An additional instability, the oscillatory skew-varicose 
(OSV) instability, is encountered for stress-free boundary 
conditions in model 1. The OSV instability boundary is 

approximately ^ = ( ~^i,^^ ) qgm, for small fi. 
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FIG. 19. (Color online) Stability diagrams in {g, q/y/JI) plane for model 2 with 7 — 2.5. (a) /i — 0.01 (b) /i = 0.001 (c) 
/i = 0.0001. In all three cases CR instability reduces the region of stable stripes which is hatched. 









Eckhaus 


0.25 


SV 


^ 
-0.25 




1 osv 


CR 


-0.5 


Existence- 




, 



0.5 



^ 



200 400 600 800 1000 

9,1, 



(a) 



Eckhaus 



0.5 



50 



100 



(b) 



Existence- 
150 200 



0.25 
> 0^ 

OH 

-0.25 : 
-0.5 



Eckhaus 







50 



100 



Existence 
150 200 



(c) 



FIG. 20. (Color online) Stability diagrams in {gm,,q/y/JJ^) plane with c — (stress-free boundary conditions), Pr — 1 and 
7 — 2.5. (a) fi = 0.1 and for large Qm- (b) fi = 0.1 and for small g-m'- stable stripes are completely eliminated when gm <; 130. 
(c) n — 0.01 and for small gm- In all three cases, the region of stable stripes (hatched ) is mainly bounded by the OSV, SV 
and CR instabilities, and the CR instability makes the upper bound in gm and reduces the region of stable stripes with fi. 



We confirmed our analytical results by numerical 
computations of the eigenvalues of the stability matri- 
ces. These eigenvalues also allow us to explore short- 
wavelength instabilities: cross-roll and the oscillatory in- 
stability. Finally, we have shown how the region of stabil- 
ity of stripes is eliminated for small fj, and large enough 
9- 

The use of the projection operator Pa, which is equiva- 
lent to a truncation to selected wave numbers, made this 
analysis straightforward and allowed the complete un- 
derstanding of the skew-varicose instability in our mod- 
els. Numerical simulations of these projected models for 
small /i have qualitatively the same solutions as the un- 
projected PDEs (for example, they have SDC solutions); 
this is our justifcation for using these projected models 
in the stability analysis for small fi. The projected and 
unprojected models will of course differ for large ^. It is 
of interest to address the question whether a similar pro- 
jection could be used in the analysis of the Navier-Stokes 
equations. 

We finally comment on the implications of our analy- 
sis on the direct numerical simulations of PDE's in large 
domains. The most striking signature of the inclusion of 



a mean-flow is the existence of the skew-varicose insta- 
bility, which can play an important role in the formation 
of the spiral defect chaos or defect chaos [l^, [iJ] ■ Hence 
the improved understanding of the stability of stripes in 
this work provides a foundation for numerical simula- 
tions of spiral defect chaos and defect chaos. The SDC 
state exists inside the Busse balloon, where convection 
rolls are stable [1[. Thus we intend in a future paper 
to relate the SDC and defect chaotic states present in 
these generalized SH models to calculations carried out 
using Rayleigh-Bcnard convection with stress-free jl| and 
no-slip [5| boundary conditions, aiming to improve the 
understanding of why SDC occurs in convection. The 
results of this work provide useful information for the 
choice of parameters for different instability regimes in 
model 1 and model 2. In particular we will be interested 
in numerical simulations with small fi and in large do- 
mains. This work also justifies using model 2 with large 
g, since this has been shown to have similar stability 
properties as to model 1 with small Pr. 
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Appendix I 

We compute the stability boundaries by defining con- 
ditions on the maxinLum eigenvalue, amax and using the 
continuation package MATCONT [3l|. 

We calculate eigenvalues of the matrices Ji and J2 and 
determine numerically how these depend on k and / in 
order to compute the derivatives that are needed. 

• Stripes are unstable to the Eckhaus instability if 
f^max is positive for / = and k = kj^ax > 0. The 
Eckhaus stability boundary occurs when kmax — ^ 0, 



in which case 



a'a„ 



afc2 



Oat (fc,0 = (0,0). 



Stripes are unstable to the zigzag instability if a 
is positive for fc = and I = Imax > 
stability boundary occurs when Imax - 



The zigzag 
0, in which 



case 



d'a- 



|p^ = Oat (fc,0 = (0,0). 



• There are two cases of the SVI, as discussed in the 
text. First if the SVI precedes the Eckhaus insta- 
bility, we consider a„iax as a function of k and / in 
polar coordinate, so A: = ecos(6') and I = esin(0). 
Stripes are unstable to the SVI, if (Tmax is positive 
for some 6 in the limit of e = 0. The SV stabil- 
ity boundary occurs when e — > 0, in which case 
a^ax = and ^^ = 0. Second if the SVI fol- 
lows the Eckhaus instability, stripes are unstable 



to the SVI, if (Jmax is positive for sufficiently small 
I = Imax and k = k^ax- The SV stability bound- 
ary occurs when Imax — > 0, in which case g'^S'' ~ 
and ^^^ = at (fc, /) = {kmax, 0). 

• In the case of CR instability, amax — occurs 
at non-zero {k,l). Hence the three conditions are 

(Pk^ ^ ' g;"^ ^ and (Tmax ^= at K ^ f^max T" 

and I = Imax 7^ 0. 



The SV and CR instabilities can be oscillatory. In 
these cases, we consider (Tmax to be the real part of the 
eigenvalue and use the same relevant conditions stated 
above. 
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